Evolutionary Trajectories of SARS-CoV-2 Alpha and Delta Variants in White-Tailed Deer in Pennsylvania (preprint)

The SARS-CoV-2 pandemic likely began by spillover from bats to humans; today multiple animal species are known to be susceptible to infection. White-tailed deer, Odocoileus virginianus are infected in the United States at substantial levels, raising concerns about the formation of a new animal reservoir and potential of spill-back of new variants into humans1. Here we characterize SARS CoV-2 in deer from Pennsylvania (PA) sampled during fall and winter 2021. Of 93 nasal swab samples analyzed by RT-qPCR, 18 (19.3%) were positive for SARS-CoV-2. Seven whole-genome sequences were obtained, which were annotated as alpha and delta variants, the first reported observations of these lineages in deer, documenting multiple new jumps from humans to deer. The alpha lineage persisted in deer after its displacement by delta in humans, and deer-derived alpha variants diverged significantly from those in humans, consistent with a distinctive evolutionary trajectory in deer.

SARS-CoV-2 variants associated with vaccine breakthrough in the Delaware Valley through summer 2021 (preprint)

The severe acute respiratory coronavirus-2 (SARS-CoV-2) is the cause of the global outbreak of COVID-19. Evidence suggests that the virus is evolving to allow efficient spread through the human population, including vaccinated individuals. Here we report a study of viral variants from surveillance of the Delaware Valley, including the city of Philadelphia, and variants infecting vaccinated subjects. We sequenced and analyzed complete viral genomes from 2621 surveillance samples from March 2020 to September 2021 and compared them to genome sequences from 159 vaccine breakthroughs. In the early spring of 2020, all detected variants were of the B.1 and closely related lineages. A mixture of lineages followed, notably including B.1.243 followed by B.1.1.7 (alpha), with other lineages present at lower levels. Later isolations were dominated by B.1.617.2 (delta) and other delta lineages; delta was the exclusive variant present by the last time sampled. To investigate whether any variants appeared preferentially in vaccine breakthroughs, we devised a model based on Bayesian autoregressive moving average logistic multinomial regression to allow rigorous comparison. This revealed that B.1.617.2 (delta) showed three-fold enrichment in vaccine breakthrough cases (odds ratio of 3; 95% credible interval 0.89-11). Viral point substitutions could also be associated with vaccine breakthroughs, notably the N501Y substitution found in the alpha, beta and gamma variants (odds ratio 2.04; 95% credible interval of 1.25-3.18). This study thus provides a detailed picture of viral evolution in the Delaware Valley and a geographically matched analysis of vaccine breakthroughs; it also introduces a rigorous statistical approach to interrogating enrichment of viral variants.

Comparative Analysis of Emerging B.1.1.7+E484K SARS-CoV-2 isolates from Pennsylvania (preprint)

Rapid whole genome sequencing of SARS-CoV-2 has presented the ability to detect new emerging variants of concern in near real time. Here we report the genome of a virus isolated in Pennsylvania in March 2021 that was identified as lineage B.1.1.7 (VOC-202012/01) that also harbors the E484K spike mutation, which has been shown to promote "escape" from neutralizing antibodies in vitro. We compare this sequence to the only 5 other B.1.1.7+E484K genomes from Pennsylvania, all of which were isolated in mid March. Beginning in February 2021, only a small number (n=60) of isolates with this profile have been detected in the US, and only a total of 253 have been reported globally (first in the UK in December 2020). Comparative genomics of all currently available high coverage B.1.1.7+E484K genomes (n=235) available on GISAID suggested the existence of 7 distinct groups or clonal complexes (CC; as defined by GNUVID) bearing the E484K mutation raising the possibility of 7 independent acquisitions of the E484K spike mutation in each background. Phylogenetic analysis suggested the presence of at least 3 distinct clades of B.1.1.7+E484K circulating in the US, with the Pennsylvanian isolates belonging to two distinct clades. Increased genomic surveillance will be crucial for detection of emerging variants of concern that can escape natural and vaccine induced immunity.

Signatures of COVID-19 severity and immune response in the respiratory tract microbiome (preprint)

Rationale: Viral infection of the respiratory tract can be associated with propagating effects on the airway microbiome, and microbiome dysbiosis may influence viral disease. Objective: To define the respiratory tract microbiome in COVID-19 and relationship disease severity, systemic immunologic features, and outcomes. Methods and Measurements: We examined 507 oropharyngeal, nasopharyngeal and endotracheal samples from 83 hospitalized COVID-19 patients, along with non-COVID patients and healthy controls. Bacterial communities were interrogated using 16S rRNA gene sequencing, commensal DNA viruses Anelloviridae and Redondoviridae were quantified by qPCR, and immune features were characterized by lymphocyte/neutrophil (L/N) ratios and deep immune profiling of peripheral blood mononuclear cells (PBMC). Main Results: COVID-19 patients had upper respiratory microbiome dysbiosis, and greater change over time than critically ill patients without COVID-19. Diversity at the first time point correlated inversely with disease severity during hospitalization, and microbiome composition was associated with L/N ratios and PBMC profiles in blood. Intubated patients showed patient-specific and dynamic lung microbiome communities, with prominence of Staphylococcus. Anelloviridae and Redondoviridae showed more frequent colonization and higher titers in severe disease. Machine learning analysis demonstrated that integrated features of the microbiome at early sampling points had high power to discriminate ultimate level of COVID-19 severity. Conclusions: The respiratory tract microbiome and commensal virome are disturbed in COVID-19, correlate with systemic immune parameters, and early microbiome features discriminate disease severity. Future studies should address clinical consequences of airway dysbiosis in COVID-19, possible use as biomarkers, and role of bacterial and viral taxa identified here in COVID-19 pathogenesis.